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Abstract. We study the out-of-equilibrium time evolution after a local quench connecting 
two anisotropic spin-1/2 XXZ Heisenberg open chains via an impurity bond. The dynamics 
is obtained by means of the adaptive time-dependent density-matrix renormalization group. 
We show that the entanglement entropies (Von Neumann and Renyi), in the presence of a 
weakened bond depend on the sign of the bulk interaction. For attractive interaction (A < 0), 
the defect turns out to be irrelevant and the evolution is asymptotically equivalent to the one 
without defect obtained by conformal field theory. For repulsive interaction (A > 0), the 
defect is relevant and the entanglement saturates to a finite value. This out-of-equilibrium 
behavior generalizes the well known results for the ground-state entanglement entropy of the 
model. 



1. Introduction 

Recent times witnessed a renewed interest in the non-equilibrium dynamics of isolated many- 
body quantum systems. A particular class of these non-equilibrium problems which is 
experiencing a dramatic explosion of theoretical activity is that of a sudden quench of a 
Hamiltonian parameter, principally boosted by the experiments on trapped ultra-cold atomic 
gases [1, 2, 3, 4, 5, 6, 7, 8] in which it has been shown that it is possible to follow 
the unitary non-equilibrium evolution without any significant coupling to the environment. 
These experiments fall into two main classes which are usually denoted as global and local 
quantum quenches. In the former case, a control parameter is suddenly quenched in the 
whole system (usually in a translationally invariant manner), while in the latter the change 
is only local, see e.g. [9] for a review. Global quantum quenches are ideal experiments 
to investigate the intriguing issue of the existence of a stationary state and thermalization 
[10, 11, 12, 13, 14, 15, 16, 17, 18], while local quenches usually reveal the spreading of 
information and correlations in a cleaner way because they mainly probe universal low- 
lying excitations. Other interesting effects can be uncovered in some 'intermediate' non- 
equilibrium situations, such as inhomogeneous quenches [19, 20, 21, 22, 23] and gas 
expansion [7, 8, 24, 25, 26, 27]. 

In this paper we consider the most commonly studied local quench in which two halves 
of a critical one-dimensional system are initially prepared in their respective ground states 
and at a given time, let us say t = 0, they are connected and let evolve according to a unitary 
Hamiltonian dynamics. Previous studies [28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38] showed 
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that a very useful tool to understand the non-equilibrium quench dynamics is the entanglement 
entropy between two complementary parts. For a general bipartition of a pure state |^) of 
quantum system (i.e. writing the whole Hilbert space of the system as a direct product of two 
parts "H = Ha ® "Hs), the Renyi entropy of the reduced density matrix pA = TrB\^){'^\ of 
the subsystem A [39] 

4") = ^lnIYp^, (1) 

is a proper measurement of the entanglement between the two parts. In the limit a — )■ 1, 
S^2^ reduces the more studied von Neumann entanglement entropy but, the knowledge of the 
Renyi entropies for any a gives far more information than the a = 1 case because it provides 
the full spectrum of the reduced density matrix [40]. 

In the ground-state of a one-dimensional critical system, whose scaling limit is described 
by a conformal field theory (CFT), in the case when A is an interval of length £ embedded 
in an infinite system, the asymptotic large i behavior of the Renyi entropies is given by 
[41,42, 43,44] 

Si'-) = ^ (l + ^]\ni + c'^, (2) 



6 \ 

where c is the central charge [45] and a non-universal constant. In the case when the whole 
system is finite, the Renyi entropies can be obtained by a conformal mapping from the cylinder 
to the plane and the net effect is just to replace i with the chordal length L/tt sin(7r£/L) [42]. 
More complicated and model dependent universal expressions for S^°'^ have been obtained for 
the low-lying excited states conformal theories [46, 47] as well. 

Conformal invariance can also be exploited to predict the behavior of the entanglement 
entropy for a 'cut and glue' local quench when two semi-infinite critical systems are joined 
together in a translational-invariant/homogeneous way (i.e. in a microscopical model, the 
added coupling is equal to all the others in the system). It has been found [29] that the 
entanglement grows logarithmically with time, again with a prefactor given by the central 
charge, i.e. 

= J: (l + - ) Int + cst. (3) 



6 \ 

This has been generalized by Stephan and Dubail [36] to the case when the two parts which 
are connected at t = hawe finite length. In the special case when the two parts have the same 
length L/2, the final result is 



L ( TlVt 

sin 



, , +cst., (4) 

where we explicitly introduced the speed of the sound v. These results, as well as some for 
general observables, have been largely tested in the literature both for free and interacting 
models [32, 33, 34, 35, 36, 48, 49, 50, 51, 52]. Furthermore it is worth mentioning that a 
few proposals to measure the entanglement entropy in real experiments are based on local 
quantum quenches [53, 54, 55, 56]. 



Entanglement evolution across defects in critical anisotropic Heisenberg chains 



3 



It is natural to wonder how the previous equilibrium and non-equilibrium results 
generalize to the case when the two halves are joined with a coupling that is not equal to 
the others, especially in view of their experimental realizations, where the control of the 
added bond can be not perfect. The simpler case of non-interacting fermions turned out to 
be very peculiar and not general enough, as opposite to the homogenous case. Indeed, it has 
been found that the defect is marginal [57, 58]. For the ground-state entanglement entropy 
this implies that the entanglement between the two parts separated by the defect still grows 
logarithmically with the system size L, but with a pre-factor Ca that continuously interpolates 
as function of the defect strength between c/6(l + 1/a) in the absence of the defect and 
when the defect is so strong to divide the system in two parts. This effect was firstly found 
numerically [59, 34] and only later the interpolating function has been exactly calculated as 
function of the defect strength [60, 61, 62], which for general a reads 



Ca{s) = ^r- r / dxln 

TT^l — a 



oo 



1 _|_ Q-2aui{x,s) 



ijj{x, s) = acosh 



cosh(x) 



(5) 



1 _^ Q-2ui{x,s)y 

where s G [0, 1] is the transmission amplitude simply related to the defect strength (see 
below). Finally, very recently, it has been shown that the same interpolating function enters 
also in the prefactor in Eq. (4) for the time evolution of the entanglement entropy after a 
local quench connecting two chains of free spinless fermions [63, 64]. All these results for 
non-interacting fermions in the presence of a defect still await a full CFT derivation which, 
also in view of the similar treatment in Ref. [66], must be clearly possible. 

For bulk interacting spinless fermions, the asymptotic behavior of the entanglement 
entropy is completely different. Indeed, in equilibrium an old renormalization group (RG) 
analysis for one-dimensional interacting fermions (modeled as a Luttinger liquid field theory) 
shows that the relevance of the defect depends on the sign of the bulk interaction [57, 58]. For 
attractive interaction, the defect turned out to be irrelevant, i.e. for large enough distances, 
the system behaves like if the defect was not present, while for repulsive interaction the 
defect is relevant and even a small defect cuts the chain into two almost independent halves. 
Also the logarithmic behavior of the entanglement entropy was calculated in Ref. [67] to the 
lowest order in the impurity strength. Finally, a numerical study of the spin-1/2 anisotropic 
Heisenberg XXZ chain with a central defect was presented in Ref. [68] and it has been 
found that the von Neumann entanglement entropy grows logarithmically with L for A < 
(ferromagnetic region corresponding to attractive interaction) with a prefactor equal to 1/3 
(i.e. confirm the irrelevance of the defect because c = 1) and instead saturates for A > 
(antiferromagnetic region corresponding to repulsive interaction), confirming the relevance 
of the impurity. The approach to the thermodynamic limit is non-universal and obviously 
depends on the defect strength. The same qualitative results have been also found for a spin-1 
gapless model [69]. 

In this paper we study the time evolution of the anisotropic XXZ Heisenberg chain after 
a local quench connecting two initially separated half-chains with a defect bond, as shown in 
Fig. 1. We tackle the problem considering the evolution of an open chain with one modified 
bond in the center and using the adaptive time-dependent density matrix renormalization 
group (tDMRG) [70]. This numerical method is well suited for the calculation of the 
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Figure 1. Graphical representation of the local quench considered in this paper: a spin-chain 
is prepared in the ground-state of two disconnected halves that at time t = are joined through 
a defect bond (represented as a spring different from the others). 



entanglement entropy since at each step of the algorithm one needs to reconstruct the reduced 
density matrix and its spectrum. The manuscript is organized as follows. In Sec. 2 we present 
the details of the model and of the numerical method. In Sec. 3 we report the time evolution 
of the entanglement (von Neuman and Renyi) entropy. The simulations have been performed 
for different values of the defect strength 6 and anisotropic parameter A, as well as for a few 
chain lengths L. Finally in the last section we draw our conclusions. 



2. The model, the local quench, and the method 

We consider the anisotropic XXZ Heisenberg chain with L sites and a defect 5 in the center 
defined by the Hamiltonian 

L-l 

H = Y,J^ [S^S:^, + S^Sl, + AStSU;\ , (6) 

2=1 

with 

and we impose open boundary conditions on the sites i = 1 and i = L. Here Sf are the 
local spin-1/2 operators, i.e. in terms of Pauli matrices = crf/2. In the homogenous 
case with 5 = 0, the model is integrable by means of Bethe ansatz [71] allowing for an exact 
characterization of the ground state and all excited states. For |A| < 1 the model is gapless, 
conformal (for A 7^ —1) and described by a Luttinger liquid field theory with central charge 
0=1 and Luttinger parameter K = 7r/(2 arccos A), while for |A| > 1 it acquires a gap. 
For A 7^ 0, although the model is Bethe ansatz solvable, it is still not known how to use 
integrability to calculate effectively the entanglement entropies in the gapless phase, in spite 
of several attempts in the literature [72] (in the gapped phase instead some exact results are 
known [42, 73]). 

Via a Jordan-Wigner transformation, apart from boundary terms, the chain in Eq. (6) is 
mapped into a lattice model of spinless fermions cj (which satisfy canonical anti-commutation 
relations {q, cl^} = 6i^rn) with Hamilltonian 

L-l 



■ =1 



1\ / + 1 
2 



(8) 
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which has attractive nearest- neighbor interaction for A < and repulsive for A > 0. The 
case A = corresponds to free spinless fermions (or XX chain). 

In this manuscript we consider the non-equilibrium situation in which the chain is 
initially prepared in the relative ground states of two equal parts of length L/2, i.e. in the 
global ground-state of the Hamiltonian (6,7) with 5o = —1. We fix L to be an integer multiple 
of 4, so that the ground state of each of the two initially separated parts is non-degenerate and 
has zero magnetization (i.e. the corresponding fermion lattice is half-filled). At time t = 0, 
these two parts are connected by a defect bond J1/2 = Jdef = 1 + 5 characterized by a deficit 
S, while the rest of the chain is left unchanged with J = 1. The transmission amplitude s in 
Eq. (5) is related to the defect strength as s = sin(2 arctan(l + 6)) [60], for —1 < 6 < 0. 
The above assumptions are not only technical: away from half-filling a marginal operator 
(corresponding to forward scattering in fermion language) can be present and it is expected 
to modify the behavior of the entanglement entropies as it is known in equilibrium for other 
observables [74, 75]. However, a proper analysis of these cases is beyond the goal of this 
paper. 

Using the general CFT result reported in the introduction, for a homogeneous quench 
(i.e. with 5 = 0) the von Neumann entanglement entropy (hereafter we will refer to S^^^ just 
as S to enlighten the notation) is [36] 



and also the anisotropy strength dependence of the sound velocity is exactly known [71] 



which is valid in the gapless phase A G [—1,1]. Notice that Eq. (9) is periodic in time with 
period tr = L/v. We will refer to tr as revival time. The existence of the revival time is a 
manifestation of the well-known light-cone effect [28, 10, 29, 36]. Indeed U is the time needed 
for an excitation moving with speed v starting from the center to arrive to the boundary, being 
reflected elastically and return to the center. In a local quench and in the thermodynamic 
limit, since the excess of energy of the initial state (compared to the ground state) is intensive, 
only low-lying universal excitations can be populated and these have all the same speed v. 
Clearly, in systems with a finite number of degrees of freedom, quasi-particles with different 
velocities could be excited even by a small energy excess resulting in imperfect revivals (the 
precise interference due to finite size-scaling effects in conformal systems is largely discussed 
in Ref. [36] to which we remand for concrete examples). Notice that since v in Eq. (10) goes 
to zero approaching A = — 1, then becomes very large and no revival will be numerically 
accessible in this limit. However, for A = — 1 the model is not conformal invariant anymore 
because the dynamical critical exponent becomes z = 2. 

2.1. Method 




(9) 



V = 



IT sin(arccos A) 



(10) 



2 arccos A 



We study the local quench dynamics of the XXZ chain by means of extensive tDMRG 
simulations [70]. The algorithm initially performs a static subroutine [76] which selects the 
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Figure 2. Comparison of the entanglement entropy between tDMRG (symbols) and 
exact results (lines) in the XX chain (A = 0) for three different impurity strengths. The 
data in the main panel agree perfectly. The inset shows the relative difference \AS\ = 
I'S'tDMRG — Scx\/Scx- The relative error is smaller than 0.2% for almost the entire duration 
of the temporal evolution. The final increase is mainly due to the decreasing of the absolute 
value of the entropy which therefore implies an increase of the relative error Nevertheless, the 
relative error is always less than 1%. 

initial state as the tensor product of the non-degenerate ground states of the two halves. In 
the decimation process, we keep a number of states such that the energy precision is at least 
of the order of 10^^. Subsequently, we perform the evolution in the presence of the impurity 
bond using the time-adapting block-decimation procedure implemented in tDMRG code [70]. 
We use the second order Suzuki- Trotter decomposition of the evolution operator with time 
step dt = 5 ■ 10~^. For each time step, the local evolution operator is applied sequentially 
on each bond starting form the left boundary of the chain. We adapt in time the number 
of states used to describe the reduced Hilbert space retaining at each local step all those 
eigenvectors of the reduced density matrix corresponding to eigenvalues larger than 10^^^, up 
to a maximum value xmax = 200 (clearly the effective maximal value used by the algorithm 
strongly depends on the simulation parameters). As well known (se e.g. Refs. [43, 77]), the 
computational complexity of the time evolution of a quantum system on a classical computer 
using the tDMRG algorithm is essentially set by the growth of the bipartite entanglement. 
As the entanglement increases with time, we have to enlarge the dimension x of the reduced 
Hilbert space in order to optimally control the truncation error. In spite of this refined adaptive 
choice of x, the truncation procedure is the main source of error of the algorithm for the largest 
system sizes, but if we would let x to grow without restrictions, the algorithm would easily 
get stuck in never ending computations. In the following we will only consider the case of a 



Entanglement evolution across defects in critical anisotropic Heisenberg chains 



7 




Figure 3. Time evolution of the entanglement entropy after connecting two XXZ spin chains 
with total length L — 64. Left: We report the results for fixed A = —0.5 and varying the 
defect strength from to —0.9 at steps of 0.1. From top to bottom, the absolute strength of 
the defect \d\ increases. Right: In this panel the defect strength is fixed to 5 = —0.4, while 
we vary the bulk anisotropy from A ~ —0.6 (topmost curve) to A = 1 (bottommost curve) at 
steps of 0.1. The dashed hne corresponds to the XX chain A = 0. 



weakened bonds, i.e. 6 G [—1, 0] in Eq. (7). 

To test the accuracy of the algorithm, we benchmarked it in the XX case (A = 0) which 
is a simple free-fermion hopping model and wherein one can exactly calculate the finite-size 
entropy from the one-particle correlation matrix {c\cj) [78, 43], even for a chain with defects 
[59]. In Fig. 2 we report the numerical results for the time evolution of the entanglement 
entropy in a XX chain of L = 64 spins after a local quench with a defect bond. In the 
figure, we confront the tDMRG data with the exact numerical calculations for three different 
values of the defect strength 5. The relative error increases with time as expected, but it 
remains bounded below 1% (in this case, this error is essentially due to the accumulation of 
the truncation errors). 

3. Time evolution of the entanglement entropy 

By means of tDMRG, we calculate numerically the time evolution of the entanglement 
entropies (von Neumann and Renyi) for spin-chains of length L = 32, 64, 128. For the smaller 
system sizes L = 32, 64 we consider many possible values of the defect strength, namely all 
values of 5 = —0.1m with m = 0, 1, . . . , 9 and 5 = —0.05, and several values of the bulk 
anisotropy parameter A G (—1, 1], namely all A = 0.1m with m = —9, —8 ... 10 (i.e. only 
values of A such that the scaling limit of the chain is described by a CFT with c = 1). 

The data for the von Neumann entanglement entropy in spin chains of length L = 64 as 
function of time are shown in Fig. 3, where in the left panel we fix A = —0.5 and report all 
the values of 5, while in the right panel we fix 5 = —0.4 and show all the considered values 
of A G [-0.6, 1] up to t = 2L/v. 
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Figure 4. Differences between the entanglement entropy evolution and the CFT prediction 
(9) in the attractive regime (A < 0) for two different values of A and four representative 
values of the defect bond. Dotted hnes are for L = 32, dashed lines for L — 64, and full lines 
for L = 128. Increasing the system size these differences tend to flatten, showing that, for 
A < 0, a defect 6 £ (—1, 0) is always irrelevant. 



For L = 128, the simulations are computationally more demanding. Thus we focus on 
few most significative values of A and 5 and we only study the non-equilibrium dynamics 
up to the revival time U = L/v since the universal part of the evolution should be periodic 
with period tr- The simulations allow us to draw the following general scenario which we 
anticipate before the careful analysis of the data. The non-equilibrium behavior after the 
local quench is reminiscent of the equilibrium one: the entropies evolve in time according 
to the CFT prediction (9), provided that A < 0, i.e. attractive interaction; conversely, for 
repulsive interactions (A > 0), the entanglement growth is suppressed and the suppression 
increases with repulsive coupling constant A. Thus, as for the equilibrium counterpart, 
there is a fundamental difference between repulsive and attractive interactions: even out- 
of-equilibrium, the operator associated to the impurity turns out to be irrelevant for A < 
and relevant for A > 0. Obviously, the 'velocity of the RG flow' toward its asymptotic limit 
strongly depends on the defect strength. In the following we report and discuss the results for 
Von Neumann and Renyi entropies to support this scenario in two separate subsections. 

3.1. Von Neumann entropy 

We now report the numerical data for the Von Neumann entanglement entropy discussing 
separately the cases of attractive and repulsive interactions. Fig. 4 shows the results for 
attractive interaction A < (note that for A = —0.8 the data for L = 128 do not reach 
the revival time L/v because this would have been computationally too demanding since v 
is very small cf. Eq. (10)). In order to show accurately that, in the thermodynamic limit, 
independently of the defect strength S E (—1,0), the entanglement entropy is described by 
the CFT prediction (9), we have subtracted the expected conformal formula (9) to the time 
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evolution of the entanglement entropy. For small (6 = —0.05) and moderate (5 = —0.4) 
defects, the difference S — Soft clearly tends to a constant value (independent on t and L) 
as long as the rescaled variable vt/L is sufficiently far from the boundaries vt/L ~ or 1, 
where finite size effects are enhanced. However, when the defect gets stronger (e.g. 6 = —0.8) 
deviations from constant behavior are evident. Increasing the system sizes, these curves tend 
to flatten for both values of the anisotropic parameter reported in the figure, suggesting that 
in the thermodynamic limit the difference would be perfectly constant for any A. A proper 
finite-size scaling analysis is made difficult by the non-uniform approach as function of vt/L, 
but already a qualitative look of the data leaves no doubts on the asymptotic result. Notice that 
for A = —0.5 larger system sizes are required for the differences to get constant compared to 
A = —0.8. Fig. 3 shows that this is a general trend as function of A, i.e. smaller is | A| larger 
is the system size required to observe the conformal behavior. The origin of this dependence 
on A is easily understood. For A = 0, the time evolution of the entanglement entropy is not 
given by the CFT formula (9), but the logarithm has a defect dependent prefactor [63]. Thus, 
for small A, increasing L there is a crossover between the A = behavior and the asymptotic 
result. 

In the case of repulsive interactions (XXZ antiferromagnet with A > 0), let us first 
discuss what we would expect if, in analogy with equilibrium, the defect would be relevant. 
For any A, the entanglement entropy should initially grow, but at some given time t* the 
presence of the relevant defect prevents further passage of information between the two halves 
and the entanglement saturates to a constant value (up to a possible revival at L/v). Clearly 
t* is not universal and its value would depend both on 5 and L. When S is small enough, 
we expect this saturation time to be so large that for any accessible size L, it would be very 
difficult to find a perceptible difference with the CFT behavior at S = 0. Increasing S, the 
saturation time t* decreases, even if for small and moderate values of L, the difference with 
the CFT evolution should remain small. Further increasing 5, the saturation time t* should 
become very short and the entanglement quickly saturates to its asymptotic value. 

If this scenario is correct, it does not make sense to subtract the conformal formula to 
tDMRG data, and so we simply report the entanglement entropy evolution for two values 
of A in Fig. 5, supplementing the data in the left panel of Fig. 3 for L = 64. From 
these two figures, it is clear that the previously conjectured scenario for the time evolution 
is confirmed by the numerics. For all values of A, the growth of the entanglement entropy 
observed in homogeneous systems is more and more suppressed as the absolute value of 
the defect strength increases. For the smallest impurity we considered (5 = —0.05), the 
entanglement entropy is very close to the conformal prediction, but the lower quality of the 
fit (that we do not report in the plot) compared to the data at 5 = suggests that some tiny 
differences can be already present. As the defect becomes stronger (e.g. S = —0.4), the curves 
flatten substantially and slowly approach an horizontal line in the thermodynamic limit. For 
very strong defect {5 = —0.8) the entanglement entropy becomes flat in a very short time. 
Of course, also in this case the bulk interaction plays a role because the saturation limit is 
approached faster for larger values of A. As in the case of attractive interactions, this is due 
to the fact that for A = there is a logarithmic growth of the entanglement for any 5 and so 
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Figure 5. Entanglement entropy evolution in the repulsive regime: the entanglement initially 
grows and then saturates to a time independent value, but the saturation time is very large for 
small 6 and would require very large systems to be observed. Dashed lines are for L = 64 and 
full lines for L = 128. 



there is a crossover in L from the non-interacting behavior to the saturation in the interacting 
chain. 

Finally we also discuss the revival of the entanglement entropy. In the right panel of Fig. 
3, we show the time evolution up to 2L/v (i.e. for two full asymptotic periods) for L = 64, 
5 = —0.4 and several values of A. For attractive interaction A < 0, the quasi-periodic 
behavior is evident, but during the second period some oscillations in time are superimposed 
to the CFT prediction (9). These have been already observed for A = (see e.g. [36]), even 
in the absence of the defect, and they are mainly due the interference between excitations with 
non perfectly equal velocities. It has been shown [36] that the amplitude of these oscillations 
goes to zero increasing the system size for A = and 5 = 0. From the figure it is evident 
that for L = 64, the amplitude of the oscillations gets smaller decreasing A, thus we expect 
they should vanish in the thermodynamic limit for all A < 0, as also corroborated by the data 
for L = 32 not shown here. For A > 0, the effect of the slower quasi-particles is different. 
Indeed, Fig. 3 shows that the interference effect gives rise to oscillations which make the 
curve flatter in time during the second period. Thus the second period entanglement entropy 
is closer to the thermodynamic limit compared to the first period value. Notice in fact that in 
Fig. 3 for A > and fixed 6 = —0.4, the average over the second period is slightly larger 
than the average over the first one, similarly as in Fig. 5 the entanglement for L = 128 is 
larger that the one for L = 64 at 5 = —0.4. 

3.2. Renyi entropies 

In this section we check that the general scenario drawn above on the basis of the von 
Neumann entropy carries over to the time evolution of the Renyi entropies for general a, 
hmiting to report data for a = 2, 3, oo (the case a = oo corresponds to the logarithm of the 
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Figure 6. Differences between the Renyi entropy evolution and the CFT prediction in the 
attractive regime for three different values of a. Dashed lines are for L — 64 and full lines for 
L = 128. We do not report the data for L = 32 to make the plots more readable. 



largest eigenvalues of the reduced density matrix, i.e. the so-called single copy entanglement 
[79]). 

In Fig. 6 we report the results for the attractive case A < 0. In analogy with the von 
Neumann entropy, we plot the difference 5*^°^ — S^^j,. For all values of a. A, and S, the 
irrelevance of the defect is confirmed by these data, indeed for small \S\ the CFT prediction 
is already approached even for relatively small values of L, while for large \S\ the difference 
tends to a constant only in the limit of large L. However, there is also an unexpected effect. 
Indeed, in Fig. 6, for small value of 5, there are subleading oscillations on top of the leading 
CFT behavior. The presence of these oscillations is well-known and it has been observed in 
the absence of the defect in all previous studies (see e.g. [32, 33, 36]) and it is explained as the 
non-equilibrium counterpart of the equilibrium unusual corrections to the scaling [80, 81] of 
the form L"^/", with K the Luttinger parameter. The frequency of these oscillations is related 
to kp (i.e. the filling) and does not dependent on the coupling strength K [80]. The interesting 
new effect is that in the presence of the defect such corrections are largely suppressed and 
almost disappear already for 6 = —0.4. This means that the CFT predictions becomes 




Figure 7. The evolution of the Renyi entropies in the repulsive regime for different value of 
a. Dashed lines are for L = QA and full lines for L — 128. 



more accurate in the presence of a defect (that however must be not too large to observe 
the flattening of the curves in Fig. 6). We stress that these oscillations are not directly related 
to the ones discussed in the previous section due to quasi-particles with velocity different from 
V, which in Fig. 6 starts manifesting only close to the revival time. 

The data for the repulsive case are reported in Fig. 7. Even in this case, the picture drawn 
from the study of the von Neumann entropy straightforwardly carries over to the Renyi: for 
the small impurities, the Renyi entropies stay very close to the conformal prediction, but as 
the defect becomes stronger, the entropies flatten and approach quickly a saturation value for 
any a. Even in the repulsive case, it is evident that the oscillating corrections to the scaling 
present for small \5\ are largely suppressed as the defect strength increases. 



4. Conclusions 



In this manuscript, we have studied the time evolution of the entanglement entropies after 
a local quench connecting two anisotropic spin- 1/2 XXZ Heisenberg chains via an impurity 
bond. By means of tDMRG, we showed that, at half-filling, the asymptotic behavior of the 
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entanglement entropies for large systems depends only on the sign of bulk interactions. For 
attractive interaction (A < 0), it turns out that the defect is irrelevant and the evolution 
is asymptotically equivalent to the conformal one without the impurity. For repulsive bulk 
interaction (A > 0), the defect becomes relevant and the entanglement entropy saturates to a 
finite value after a finite time which are both not universal and depend on the defect strength. 

We have only determined time evolution of the entanglement entropy between the two 
initially disconnected parts, but we can safely conclude that the same picture should carry over 
to all other observables which can be determined just by exploiting the relevance/irrelevance 
of the defect according to the sign of the bulk interaction. We mention that the correspondence 
between the equilibrium and out of equilibrium situations resulting from a local perturbation 
for observables different from the entanglement entropies has a long history [74, 75] such 
as the connection between equilibrium Anderson orthogonality catastrophe and the X-ray 
absorption spectra (corresponding to a local perturbation). Furthermore, while we have 
only considered the XXZ spin chain, in view of the Luttinger liquid universality, the same 
conclusions apply to all other models in the same universality class, such as one-dimensional 
gases of spinless bosons (both the continuum Lieb-Liniger and the Bose-Hubbard model in the 
superfluid phase), quantum wires junction etc. Oppositely, for multi-component continuum or 
lattice models of interacting fermions (such as the Hubbard model) or Fermi-Bose mixtures, 
the phenomenon of spin-charge separation could lead to different universal features, e.g. 
reminiscent of what observed in Ref. [82] for different initial states. For these multi- 
component models, also shell-filling effects [47] could strongly affect the non-equilibrium 
dynamics of the entanglement entropy. Therefore it would be really interesting to investigate 
local quenches in these models both in the absence and in the presence of defect bonds. 

Acknowledgments 

Both authors acknowledge the financial support by the ERC under Starting Grant 279391 
EDEQS. 

References 

[1] M. Greiner, O. Mandel, T. W. Hansch, and I. Bloch, Nature 419 51 (2002). 
[2] T. Kinoshita, T. Wenger, D. S. Weiss, Nature 440, 900 (2006). 

[3] S. Hofferberth, I. Lesanovsky, B. Fischer, T. Schumm, and J. Schmiedmayer, Nature 449, 324 (2007). 
[4] S. Trotzky Y.-A. Chen, A. Flesch, I. P. McCulloch, U. Schollwock, J. Eisert, and I. Bloch, Nature Phys. 8, 
325 (2012). 

[5] M. Cheneau, P. Barmettler, D. Poletti, M. Endres, P. Schauss, T. Fukuhara, C. Gross, I. Bloch, C. Kollath, 

and S. Kuhr, Nature 481, 484 (2012). 
[6] M. Gring, M. Kuhnert, T. Langen, T. Kitagawa, B. Rauer, M. Schreitl, I. Mazets, D. A. Smith, E. Demler, 

and J. Schmiedmayer, Science 337, 1318 (2012). 
[7] U. Schneider, L. Hackermiiller, J. P. Ronzheimer, S. Will, S. Braun, T. Best, I. Bloch, E. Demler, S. Mandt, 

D. Rasch, and A. Rosch, Nature Phys. 8, 213 (2012). 
[8] J. P. Ronzheimer, M. Schreiber, S. Braun, S. S. Hodgman, S. Langer, I. P. McCulloch, F. Heidrich-Meisner, 

I. Bloch, and U. Schneider, ai-Xiv: 1301.5329. 
[9] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, Rev. Mod. Phys. 83, 863 (2011). 



Entanglement evolution across defects in critical anisotropic Heisenberg chains 



14 



[10] P. Calabrese and J. Cardy, Phys. Rev. Lett. 96, 136801 (2006); 

P Calabrese and J. Cardy, J. Stat. Mech. (2007) P06008. 
[1 IJ M. Rigol, V. Dunjko, V. Yurovsky, and M. Olshanii, Phys. Rev. Lett. 98, 50405 (2007); 

M. Rigol, V. Dunjko, and M. Olshanii, Nature 452, 854 (2008). 
[12] M. Cramer, C. M. Dawson, J. Eisert, and T. J. Osborne, Phys. Rev. Lett. 100, 030602 (2008); 

M. Cramer and J. Eisert, New J. Phys. 12, 055020 (2010). 
[13] T. Barthel and U. SchoUwock, Phys. Rev. Lett. 100, 100601 (2008). 
[14] P Calabrese, FH.L. Essler and M. Fagotti, Phys. Rev. Lett. 106, 227203 (201 1); 

P Calabrese, EH.L. Essler and M. Fagotti, J. Stat. Mech. (2012) P07016 

P Calabrese, EH.L. Essler and M. Fagotti, J. Stat. Mech. (2012) P07022. 
[15] R H. L. Essler, S. EvangeUsti, and M. Fagotti, Phys. Rev. Lett. 109, 247206 (2012). 
[16] D. Schuricht and F H. L. Essler, J. Stat. Mech. (2012) P04017. 
[17] G. Biroh, C. Kollath, and A. Laeuchli, Phys. Rev. Lett. 105, 250401 (2010); 

G. P Brandino, A. De Luca, R.M. Konik, and G. Mussardo, Phys. Rev. B 85, 214435 (2012). 
[18] D. Fioretto and G. Mussardo, New J. Phys. 12, 055015 (2010); 

S. Sotiriadis, D. Fioretto, and G. Mussardo, J. Stat. Mech. (2012) P02017. 
[19] P Calabrese, C. Hagendorf, and P Le Doussal, J. Stat. Mech. P07013 (2008). 
[20] S. Sotiriadis and J. Cardy, J. Stat. Mech. PI 1003 (2008). 
[21] V. Eisler, F Igloi, and I. Peschel, J. Stat. Mech. P02011 (2009). 
[22] T. Antal, Z. Racz, A. Rakos, and G. M. Schutz, Phys. Rev. E 59, 4912 (1999); 

D. Karevski, Eur. Phys. J. B 27, 147 (2001); 

Y. Ogata, Phys. Rev. E 66, 066123 (2002); 

T. Platini and D. Karevski, Eur Phys. J. B 48, 225 (2005); 

T. Platini and D. Karevski, J. Phys. A 40, 1711 (2007); 

T. Antal, P L. Krapivsky, and A. Rakos, Phys. Rev. E 78, 061115 (2008); 

J. Lancaster and A. Mitra, Phys. Rev. E 81, 061 134 (2010); 

S. Langer, M. Heyl, I. P McCulloch, and R Heidrich-Meisner, Phys. Rev. B 84, 205 1 15 (201 1); 

M. CoUura, H. Aufderheide, G. Roux, and D. Karevski; Phys. Rev. A 86, 013615 (2012); 

V. Eisler and Z. Racz, arXiv: 121 1.2321. 
[23] J. Mossel, G. Palacios, and J.-S. Caux, J. Stat. Mech. L09001 (2010). 
[24] A. Minguzzi and D.M. Gangardt, Phys. Rev. Lett. 94, 240404 (2005); 

D. M. Gangardt and M. Pustilnik, Phys. Rev. A 77, 041604 (2008). 
[25] M. Campostrini and E. Vicari, Phys. Rev. A 82, 063636 (2010); 

E. Vicari, Phys. Rev. A 85, 062324 (2012); 
J. Nespolo and E. Vicari, arXiv: 1301.3693. 

[26] J.-S. Caux and R. M. Konik, Phys. Rev. Lett. 109, 175301 (2012); 

G. Brandino, J.-S. Caux, and R. M. Konik, arXiv: 1301.0308. 
[27] F. Heidrich-Meisner, M. Rigol, A. Muramatsu, A.E. Feiguin, E. Dagotto, Phys. Rev. A 78, 013620 (2008); 

S. Langer, F. Heidrich-Meisner, J. Gemmer, LP McCulloch, U. SchoUwoeck, Phys. Rev. B 79, 214409 

(2009); 

G. Roux, Phys. Rev. A 81, 053604 (2010). 
[28] P Calabrese and J. Cardy, J. Stat. Mech. (2005) P04010. 
[29] P Calabrese and J. Cardy, J. Stat. Mech. P10004 (2007). 
[30] A. Laeuchli and C. Kollath, J. Stat. Mech. (2008) P05018. 
[31] M. Fagotti and P Calabrese, Phys. Rev. A 78, 010306 (2008). 
[32] V. Eisler and I. Peschel, J. Stat. Mech. (2007) P06005. 

[33] V. Eisler, D. Karevski, T. Platini, and I. Peschel, J. Stat. Mech. (2008) P01023. 

[34] F Igloi, Z. Szatmari, and Y.-C. Lin, Phys. Rev. B 80, 024405 (2009). 

[35] J.-M. Stephan and J. Dubail, J. Stat. Mech. (2011) L03002. 

[36] J.-M. Stephan and J. Dubail, J. Stat. Mech. (2011) P08019 

[37] R Igloi, Z. Szatmari, and Y.-C. Lin, Phys. Rev. B 85 094417 (2012). 



Entanglement evolution across defects in critical anisotropic Heisenberg chains 



15 



[38] G. C. Levine, M. J. Bantegui, and J. A. Burg, Phys. Rev. B 86, 174202 (2012). 
[39] L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. Phys. 80, 517 (2008); 

J. Eisert, M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 82, 277 (2010); 

R Calabrese, J. Cardy, and B. Doyon Eds, J. Phys. A 42 500301 (2009). 
[40] P Calabrese and A. Lefevre, Phys. Rev. A 78, 032329 (2008). 
[41] C. Holzhey, F. Larsen, and F. Wilczek, Nucl. Phys. B 424, 443 (1994). 
[42] R Calabrese and J. Cardy, J. Stat. Mech. P06002 (2004). 

[43] G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. Lett. 90, 227902 (2003); 

J. I. Latorre, E. Rico, and G. Vidal, Quant. Inf. Comp. 4, 048 (2004). 
[44] R Calabrese and J. Cardy, J. Phys. A 42, 504005 (2009). 
[45] J. Cardy, J. Stat. Mech. (2010) P10004. 

[46] F. C. Alcaraz, M. Ibanez Berganza, and G. Sierra, Phys. Rev. Lett. 106, 201601(2011); J. Stat. Mech. 
(2012) P01016. 

[47] F. H. Essler, A. LauchU, and P. Calabrese, Phys. Rev. Lett, to appear [arXiv:1211.2474]. 
[48] D. Gobert, C. KoUath, U. SchoIIwoeck, and G. Schuetz, Phys. Rev. E 71, 036102 (2005). 
[49] A. Silva, Phys. Rev. Lett. 101, 120603 (2008). 
[50] U. Divakaran, R Igloi, and H. Rieger, J. Stat. Mech. P10027 (2011) 

[51] M. Ganahl, E. Rabel, R H. L. Essler, and H. G. Evertz, Phys. Rev. Lett. 108, 077206 (2012). 

[52] P Smacchia and A. Silva, Phys. Rev. Lett. 109, 037202 (2012). 

[53] 1. Klich andL. Levitov, Phys. Rev. Lett. 102, 100502 (2009). 

[54] B. Hsu, E. Grosfeld, and E. Rradkin, Phys. Rev. B 80, 235412 (2009). 

[55] J. Cardy, Phys. Rev. Lett. 106, 150404 (2011). 

[56] D. A. Abanin and E. Demler, Phys. Rev. Lett. 109, 020504 (2012). 

[57] C. L. Kane and M. P Risher, Phys. Rev. Lett. 68, 1220 (1992); 

C. L. Kane and M. R Risher, Phys. Rev. B 46, 15233 (1992). 
[58] S. Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992). 
[59] I. Peschel, J. Phys. A 38, 4327 (2005). 
[60] V. Eisler and L Peschel, Ann. Phys. (Berhn) 522 (2010) 679. 
[61] R Calabrese, M. Mintchev, and E. Vicari, J. Phys. A 45 (2012) 105206; 

P Calabrese, M. Mintchev, and E. Vicari, Phys. Rev. Lett. 107, 020601 (2011). 
[62] I. Peschel and V. Eisler, J. Phys. A 45, 155301 (2012). 
[63] V. Eisler and L Peschel, EPL 99, 20001 (2012). 
[64] B. Bertini, P. Calabrese, and M. Mintchev, in preparation. 
[65] V. Eisler and S. Garmon. Phys. Rev. B 82, 174202 (2010). 
[66] K. Sakai and Y. Satoh, JHEP 0812: 001 (2008); 

C. Bachas, J. de Boer, R. Dijkgraaf, and H. Ooguri, JHEP 0206, 027 (2002). 
[67] G. C. Levine, Phys. Rev. Lett. 93, 266402 (2004). 
[68] J. Zhao, I. Peschel and X. Wang, Phys. Rev. B 73, 024417 (2006). 
[69] J. Ren, S. Zhu, and X. Hao, J. Phys. B 42 (2009) 015504. 
[70] S. R. White and A. E. Reiguin, Phys. Rev. Lett. 93, 076401(2004); 

A. J. Daley, C. Koflath, U. Schollwock, and G. Vidal, J. Stat. Mech. (2004) P04005. 

[71] V. E. Korepin, N. M. Bogoliubov and A. G. Izergin, Quantum Inverse Scattering Method and Correlation 

Functions, Cambridge University Press (1993). 
[72] J. Sato, M. Shiroishi, M. Takahashi, J. Stat. Mech. P12017 (2006); 

J. Sato and M. Shiroishi, J. Phys. A 40, 8739 (2007); 

J. L. Jacobsen and H. Saleur, Phys. Rev. Lett. 100, 087205 (2008); 

J. Damerau, R. Gohmann, N. P. Hasenclever, and A. Kliimper, J. Phys. A 40, 4439 (2007); 

B. Nienhuis, M. Campostrini, and R Calabrese, J. Stat. Mech. (2009) P02063; 
L. Banchi, R Colomo, and R Verrucchi, Phys. Rev. A 80, 022341 (2009); 

V. Alba, M. Ragotti, and R Calabrese, J. Stat. Mech. (2009) P10020; 
H. Katsura and L Maruyama, J. Phys. A 43 (2010) 175003. 



Entanglement evolution across defects in critical anisotropic Heisenberg chains 



16 



[73] R. Weston, J. Stat. Mech. L03002 (2006); 

E. Ercolessi, S. Evangelisti, and F. Ravanini, Phys. Lett. A 374, 2101 (2010); 

R Calabrese, J. Cardy, and I. Peschel, J. Stat. Mech. P09003 (2010); 

E. Ercolessi, S. Evangelisti, F. Franchini, and F. Ravanini, Phys. Rev. B 83, 012402 (2011); 

E. Ercolessi, S. Evangelisti, F. Franchini, and F. Ravanini, Phys. Rev. B 85, 115428 (2012); 

0. A. Castro- Alvaredo and B. Doyon, J. Stat. Mech. P02001 (2011). 
[74] I. Affleck and A. W. W. Ludwig, J. Phys. A 27 (1994) 5375. 

[75] A. M. Zagoskin and 1. Affleck, J. Phys. A 30 (1997) 5743. 
[76] S. R. White, Phys. Rev Lett. 69, 2863 (1992); 

S. R. White, Phys. Rev. B 48, 10345 (1993); 

U. SchoUwock, Rev. Mod. Phys. 77, 259 (2005). 
[77] J. I. Cirac and F. Verstraete, Renormalization and tensor product states in spin chains and lattices, J. Phys. 

A 42, 504004 (2009); G. Vidal, Entanglement Renormalization: an introduction, in Understanding 

Quantum Phase Transitions, ed. by L. D. Carr (Taylor & Francis, Boca Raton, 2010) arXiv:0912.1651. 
[78] I. Peschel, J. Phys. A 36, L205 (2003); 

1. Peschel, J. Stat. Mech. (2004) P06004; 

I. Peschel and V. Eisler, J. Phys. A 42, 504003 (2009); 
I. Peschel, Braz. J. Phys. 42, 267 (2012). 
[79] J. Eisert and M. Cramer, Phys. Rev. A 72, 0421 12 (2005); 
I. Peschel and J. Zhao, J. Stat. Mech. PI 1002 (2005); 

R. Orus, J.I. Latorre, J. Eisert, and M. Cramer, Phys. Rev. A 73, 060303 (2006). 
[80] R Calabrese, M. Campostrini, F. Essler, and B. Nienhuis, Phys. Rev. Lett. 104, 095701 (2010); 
J. Cardy and P Calabrese, J. Stat. Mech. (2010) P04023; 
P Calabrese and F H. L. Essler, J. Stat. Mech. (2010) P08029; 
J. C. Xavier and R C. Alcaraz, Phys. Rev. B 83, 214425 (2011); 

M. Dalmonte, E. Ercolessi, L. Taddia, Phys. Rev B 84, 0851 10 (201 1); Phys. Rev. B 85, 165112 (2012); 

P Calabrese, M. Mintchev, and E. Vicari, J. Stat. Mech. P09028 (201 1). 
[81] N. Laflorencie, E. S. Sorensen, M.-S. Chang, and I. Affleck, Phys. Rev. Lett. 96, 100603 (2006); 

M. Fagotti and P Calabrese, J. Stat. Mech. P01017 (2011). 
[82] A. Kleine, C. KoUath, L P McCuUoch, T. Giamarchi, and U. SchoUwock, New J. Phys. 10, 045025 (2008). 



